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SUMMARY 


This report examines the refraction effects of the earth’s atmosphere on the NASA 136-MHz 
Minitrack interferometer during aircraft calibration and spacecraft tracking. It is shown that only the 
troposphere (sea level to 30 km) affects aircraft calibration. It is also seen that the ionosphere (85 km 
to 1000 km) is the principal refraction bias source during spacecraft track. This results from the tro- 
posphere correction inherent in an interferometer angle calculation. 

The bias induced by the daytime ionosphere at 25° elevation can be as large as 1 mrad. This bias 
decreases to 0. 1 mrad at 80° elevation. During nighttime, when total ionosphere electron content de- 
creases, refraction generally drops to only 10% of the daytime bias. Possible exceptions have been 
cited for near-equatorial stations where anomalous ionospheric behavior during nighttime has caused 
signal dropouts of Minitrack telemetry and tracking data. When present, this anomaly has been noted 
generally to peak at local midnight. If daytime data below elevation angles of 45° are to be used for 
precision orbit computation, it is clear that ionospheric refraction corrections are a necessary part of 
Minitrack data processing. These corrections are necessary to reduce the ionospheric refraction bias to 
a value comparable to the Minitrack system resolution of 0. 1 mrad. 
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NASA MINITRACK INTERFEROMETER 


REFRACTION CORRECTIONS* 


by 

P. E. Schmid 

Goddard Space Flight Center 


1.0 INTRODUCTION 

This report reviews the influence of the earth’s atmosphere upon NASA Minitrack measurements 
during system calibration and spacecraft tracking. The need for such a review stems from the in- 
creased number of precision orbits now being based upon Minitrack data (e.g., Reference 1), coupled 
with a lack of formal documentation in the area of Minitrack refraction modeling. The basic consider- 
ations presented in this report apply equally well to any short baseline radio interferometer tracking 
system. 

Section 2.0 of this report briefly compares the angle information provided by an interferometer 
such as Minitrack with that obtained from steerable antenna systems. NASA steerable, or monopulse, 
tracking systems include— 

(1) The GSFC Range and Range Rate system (Reference 2). 

(2) The Unified S-Band system (Reference 3). 

(3) The Pulse C-Band radar (Reference 4). 

(4) The Applications Technology Satellite (ATS) Range and Range Rate system (Reference 5). 

The earth’s atmosphere can, for purposes of radiowave propagation analysis, be considered to 
consist of three layers. These layers, as indicated by Figure 1, are the troposphere (0 to 30 km), the 
free-space region (30 to 85 km), and the ionosphere (85 to 1000 km). Section 3.0 presents a brief 
description of the frequency independent tropospheric angle correction equations appropriate to steer- 
able antenna spacecraft tracking. Section 4.0 indicates why the interferometer, to a first order, inher- 
ently corrects for tropospheric refraction when the baseline length is measured in free-space wave- 
lengths at the spacecraft transmitter frequency. In section 5.0, it is pointed out that the effect of the 
one-way aircraft Doppler shift upon Minitrack calibration is negligible. For high-accuracy orbits, 
spacecraft Doppler effects (i.e., no effect at zenith, approximately 0.05-mrad maximum at £ = 30°) 
should be modeled, and easily can be, in the orbit computation program. 

Section 5.0 of this report reviews the philosophy of the present Minitrack aircraft calibration 
procedure, indicating the significance of radio and optical refraction. Equations are derived that 

*This report is a republication of NASA-GSFC Document X-55 1-69-434, October 1969. 
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TARGET ABOVE IONOSPHERE 
TRANSMITTING AT FREQUENCY f 



n= l + NdO®) 

TROPOSPHERE N(h)=Nse‘'<h 

WHERE: k = 0.16 km"' 

-40.3 Ng(h) 

IONOSPHERE N(h)= ^ — I0« 

where: Ne= ELECTRON DENSITY 
(ELECTR0NS/m3) 

f= FREQUENCY (Hz) 


Figure 1 —Relationship between troposphere and ionosphere. 


properly model the atmosphere for purposes of calibration computations. These particular equations 
are incorporated in the Minitrack star background photoplate-reduction computer program used by 
the Physical Science Laboratory, New Mexico State University, for the GSFC. 

Finally, section 6.0 indicates what improvements can be anticipated by modeling the ionosphere. 
Ray bending decreases in the ionosphere in proportion to the square of the transmission frequency 
(Reference 6). Thus, with an interferometer operating at S-Band frequencies (nominally 2 GHz), the 
distortion introduced by the ionosphere is reduced by a factor of approximately 200 relative to that at 
136 MHz. By contrast, the tropospheric bending is independent of the frequency in the range of 100 
MHz to 20 GHz (Reference 7). 

Figure 2 indicates the overall Minitrack data handling and aircraft calibration interface. With re- 
gard to Figure 2, this report is primarily concerned with the optical and radio refraction corrections 
applied to the aircraft calibration data and with the inherent tropospheric corrections associated with 
data delivered to the orbit computation user and ionospheric refraction during spacecraft tracking. 

Currently, aircraft calibration flights are performed approximately twice per year for each of the 
NASA Minitrack stations. (Site locations are presented in Figure 3.) The Minitrack and the GSFC 
Range and Range Rate systems are the primary tracking systems associated with the NASA Satellite 
Tracking and Data Acquisition Network (STADAN). 

It should be mentioned that the primary purpose of the aircraft calibration is to determine the 
far-field phase characteristics of the Minitrack antenna array. The possibility of obtaining the same in- 
formation via analysis of residuals from actual Minitrack satellite orbit calculations is currently under 
investigation. * A regression analysis of this type appears feasible if multiple station, short arc passes 


*Oosterhout, J.D., Code 514 GSFC, “Minitrack Calibration,” memorandum to C.H. Looney, Jr., Code 540 GSFC, August 19, 1968. 
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Figure 2— Overall Minitrack data handling. 
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Figure 3— Satellite Tracking and Data Acquisition Network (ST AD AN) Minitrack stations.* 


*Goddard Directory of Tracking Station Locations, NASA-GSFC Document X-554-67-54 as updated through November 1968. 
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are considered. Uncertainties in the earth geopotential function tend to mask tracking system biases 
when more than one earth revolution is calculated (Reference 8). 

2.0 INTERFEROMETER VERSUS STEERABLE 

ANTENNA ANGLE MEASUREMENT 

The Minitrack system is an interferometer tracking scheme (References 9 and 10) that provides 
angular data in the form of direction cosines describing the position of a spacecraft relative to a given 
tracking station. It is a passive system in which the spacecraft transmits at a nominal 136 MHz. Hence 
no ranging information is available, and at present no attempt is made to utilize the available, one-way 
Doppler frequency information for range rate calculations. The relationships between direction 
cosines and the conventional definitions of azimuth and elevation are given in Figure 4. As will be 
shown, the Minitrack direction cosine calculation does not provide a measure of either signal elevation 
angle of arrival or line-of-sight elevation to spacecraft but rather an elevation angle between the two, 
which, to a first order, is the angle of arrival corrected for tropospheric refraction. Elevation angle of 
arrival is the elevation angle of the incoming, essentially planar wavefront. Line-of-sight elevation is 
with reference to the straight line segment joining the interferometer array center and the spacecraft. 

This is in contrast to the elevation angle-of-arrival measure provided by tracking systems employ- 
ing steerable antennas (e.g., the NASA Unified S-Band, GSFC Range and Range Rate, and Pulse C-Band 
radar systems). The elevation angle of arrival is of course available directly from an azimuth-elevation 
antenna mount. It can be obtained by way of straightforward coordinate transformations when the 
angle data are measured with XY or polar antenna mounts (Reference 11). 

3.0 TROPOSPHERIC ANGLE CORRECTION APPLIED 

TO STEERABLE ANTENNA TRACKING SYSTEMS 

The usual correction scheme applied to elevation angle data obtained from tracking systems em- 
ploying steerable antennas is 


E = M - AE. , 


( 1 ) 


where 

E^ is the observed elevation angle of arrival, 

AE is a correction accounting for atmospheric bending, 

AE^ is associated with tropospheric bending, 
and 

AE. is associated with ionospheric bending. 

It should be pointed out that the effect of the ionosphere’s presence is always to increase the error 
angle AE above that due to the troposphere alone. Since the bending within the ionosphere above the 
height of maximum ionization (h = 300 km) is in a direction opposite to that within the troposphere, 
one might assume that the effects are partially compensatory. That this is not the case will be shown in 
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cos msCOS^S 


= SIN AZ COS E 
= COS AZ COS E 


Figure 4— Definition of Minitrack coordinate system. 


section 6.0. However, an intuitive answer may be obtained with reference to Figure 5. By considering 
a fixed angle of arrival one can observe what the error EE will be for various spacecraft altitudes. 

For a given angle of arrival, the ray path is determined; hence, the line-of-sight angle E will always 
be between the tracking station and a point on the ray trace corresponding to the spacecraft altitude /j. 
As is suggested by Figure 5, the error angle (A£ = E^ - E) will never be less than the error due to the 
troposphere alone. The error actually has a maximum when the spacecraft is within the ionosphere; it 
approaches that due to the troposphere alone as the spacecraft reaches an altitude of several thousand 
kilometers above the earth’s surface. 

The correction represented by EE^ (Equation 1) is frequency-independent in the range 100 to 
20,000 MHz. The correction represented by EE^ (Equation 1) varies inversely as the square of fre- 
quency. At frequencies greater than 1 GHz and at elevations greater than 5°, the ionosphere accounts 
for no more than 2% of the total refraction (Reference 6). At 136 MHz, however, the tropospheric 
and ionospheric refraction effects can be comparable, depending upon spacecraft altitude. At altitudes 
greater than 100 km and for >10°, the first-order tropospheric elevation correction for tracking 
spacecraft with a steerable antenna is given by (References 12, 13, and 14) 

=A/10-6)cot£„, (2) 

where and E^ are as defined previously and N^, the surface refractivity, is a convenient measure of 
surface index of refraction: 


1+A^^(10-®). 


(3) 
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Ep = ANGLE OF ARRIVAL (FIXED FOR A GIVEN RAY PATH) = CONSTANT 
E = LINE-OF-SIGHT ELEVATION (E , = LINE-OF-SIGHT ELEVATION TO SPACECRAFT AT h, 
AND Ej = LINE-OF-SIGHT ELEVATION TO SPACECRAFT AT hj) 

ERROR ANGLE = AE -- Eg ~ E 

NOTE THAT AE , > AEj AND AE-* CONSTANT AS 


Figure 5— Increase in elevation angle error due to ionosphere. 


The surface refractivity can be computed from 


77 6 / 4810c RH 

N =-^(p + ^ 

S T \ T 


(4) 


where 

T = temperature (degrees Kelvin), 

P = total atmospheric pressure (millibars), 

RH = relative humidity, 
and 

Cj = saturation vapor pressure (millibars). 

Equations 1 and 2 can also be written in terms of direction cosines. For frequencies above 1 GHz, the 
ionospheric contribution to ray bending is usually considered negligible and Equation 1 becomes 
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E = E„ - AT f 10"^) cot £ , 
o ' o ^ 

or 

cos£'= cos -A^^(10"®)cot£'J = [l +N^(10"^)] cosE^ , 


since, for E^ > 5°, 


iV/ 10-6) cot <E^. 


(5) 

( 6 ) 


However, from Figure 4, 

/ = cos a = sin AZ cos E 
m = cos j3 = cos AZ cos E . 

If and represent angle-of-arrival direction cosines and I, m represent “best estimate” of line-of- 
sight direction cosines, it is seen that: 


and 


[l+iV^(10-6)]/^ (7) 

m= [l+iV^(10-6)]m^. (8) 


In summary then, it is seen that a steerable antenna (i.e., monopulse) tracking system, when cor- 
rected for first-order tropospheric bending effects, results in the angle-of-arrival direction cosine ad- 
justments indicated by Equations 7 and 8. The next section will show that this tropospheric correction 
is inherent in the Minitrack direction cosine calculation. 

4.0 TROPOSPHERIC ANGLE CORRECTION INHERENT 

IN INTERFEROMETER MEASUREMENTS 

An ideal radio interferometer measurement will be defined as one where the incident signal is 
planar and operation is within a stable (i.e., time-invariant over measurement period), spherically strati- 
fied, exponential troposphere. These conditions are, for all practical purposes, met when radio energy 
is received from a spacecraft at altitudes in excess of 1 00 km above the earth’s surface and when fre- 
quencies in excess of 1 GHz are employed, thus obviating consideration of the ionosphere. At lower 
frequencies, such as the 136-MHz transmission of Minitrack signals, the ionosphere will also affect the 
angle measurement. However, this ionospheric effect can be shown to produce a linear addition to 
elevation angle measurement which will appear as a residual error once the troposphere has been taken 
into account. The following derivation shows that the ideal interferometer provides a built-in first 
order correction for the tropospheric effect and that this correction is the same as is normally applied 
to elevation angle data obtained from steerable antenna trackers (i.e.. Equation 2). 

4.1 The Basic Minitrack Measurement 

As suggested by Figure 6, the basic measurements of an orthogonally intersecting interferometer, 
such as Minitrack, are the direction cosines / and m. More precisely, however, the measurement is a 
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phase comparison, which, for each antenna pair, provides a measure of radio path length difference. 
This is indicated in Figure 6, where S 2 -S 1 and S^-S^ in actual operation represent radio path length 
differences rather than geometric differences associated with operation within free space. The object 
of all elevation angle corrections to radio tracking data is to approximate the geometric line-of-sight 
elevation angle between tracking stations and spacecraft as closely as possible. This line of sight is 
indicated by S in Figure 6. 

The phase measurement associated with the two interferometer antenna pairs is indicated in 
Figure 7. 

With reference to Figure 6, the direction cosines and of the angles of arrival are given by 




lo = cos 

II 

(9) 

= cos 

^2\ 

(10) 

b ’ 


where and M 2 are whole and fractional cycles of phase shift measured by the Minitrack system, 
is the signal wavelength which would be measured in the refractive media in the immediate vicinity of 
the interferometer, and b is the baseline separation of each of the interferometer legs. 

4.2 Minitrack Angle Computation 

The information contained in the raw tracking data is in terms of fractional cycles of phase shift. 
At GSFC, ambiguities are resolved and the phase counts are converted to direction cosines with baseline 
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Figure 6— Simplified diagram of Minitrack interferometer spacecraft 

tracking. 



WHERE Oq and ySgARE RESPECTIVE INTERFEROMETER 
ANGLES CORRESPONDING TO THE WAVEFRONT ANGLE 
OF ARRIVAL , E o 

Figure 7— Interferometer phase measurement. 


separation in free-space wavelengths rather than the actual wavelength indicated in Equations 9 and 
10 (Reference 1 5). Thus the angles a and |3' calculated with free-space baseline separation are not a 
measure of angle of arrival but some “pseudo angle” that is linked to the angle of arrival by the surface 
index of refraction. That is, 



and hence 

cos a' = /' = cos = [l 

Similarly, 

w'= [l +iV^(10'^)]m^ , (13) 


where /' and m' are calculated direction cosines and and are angle of arrival direction cosines. 
However, Equations 12 and 13 are seen to be exactly the results presented previously by Equations 7 
and 8 where a first-order tropospheric correction had been applied to a monopulse (steerable antenna) 
tracking system. Thus the Minitrack angle calculation has a built-in tropospheric refraction correction 
which, for £'^>10°, accounts for at least 98% of the tropospheric bending. If no atmospheric correc- 
tions are applied to Minitrack data it can be expected that the primary refraction effect of the earth’s 
atmosphere will be linked to the ionosphere. 


5.0 OPTICAL AND RADIO REFRACTION 
DURING MINITRACK CALIBRATION 

At present the Minitrack system is calibrated by means of aircraft that carry both optical and 
radio beacons. The optical beacon is photographed against a star background which permits a precise 
means for determining the line of sight between station and aircraft (Reference 16).* 

The primary purpose of the aircraft flyby is to determine the far-field phase characteristics of the 
Minitrack antenna array. Since Minitrack antenna spacing corresponding to each of the intersecting 


*“Minitrack System Training Manual,” (revised September 1958) compiled by the Tracking and Guidance Branch, Project 
Vanguard, National Aeronautics and Space Administration, Washington, D.C. Part VI, Calibration of Minitrack, by Habib, E. J., Bradford, 
J. N., Herbert, J. H., Engels, P. D., and Oosterhout, J. D. Can be obtained from the STADAN Technical Reference Control Facility, 
Code 533, Goddard Space Flight Center, Greenbelt, Maryland. 
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baselines is on the order of 50 wavelengths, the far field should be probed at slant ranges of approxi- 
mately 1 0 km or greater. The aircraft also provide a zero-set phase measurement that accounts for in- 
equalities in cable electrical length between respective antenna arrays and the Minitrack receiver. These 
cables are filled with moisture-free nitrogen and maintained at an approximate gauge pressure of 5 psi. 
In addition to the aircraft flyby, internal Minitrack system calibrations are performed to account for 
electronic equipment phase delays. 

5.1 Calibration and Optical Refraction 

As indicated in Figure 8, the calibration aircraft is photographed against a star background. The 
line of sight from the center of the Minitrack’s orthogonal pair of interferometers to a given star is 
then accurately determined since the star field has been well established by astronomers. However, the 
line of sight to the aircraft differs from the line of sight to the star by a small angle e. This angle, e, 
must be added to the line-of-sight angle to the star to provide an optical line of sight to the aircraft 
which is then corrected for radio refraction to permit a calculation of what the interferometer should 
read at 136 MHz. An expression for e has been rigorously derived in Reference 1 7. As will be shown 
in section 5.3 it is also given to a good approximation by 

e = (n - 2 + — I cot E , 

\ ^o) 

^{\(T^)NcotE (14) 


where 


E = elevation angle to the star 
and 

n = surface optical index of refraction. 

The term n , the average optical index of refraction between Minitrack station and aircraft altitudes, 
is derived from 



where 


= altitude of aircraft above sea level, 

H = altitude of tracking station above sea level, 

and 

n = optical index of refraction as a function of height, 

which for a spherically stratified, exponential troposphere, where h - ~ [^i.e., n{h) = 1 + 

7V^(10'6)e"T'^ , leads to 

(10”^) cot £■ [l - 


(15) 
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Figure 8-Optical refraction geometry. 


where the decay constant 7 for the optical index of refraction within the earth’s atmosphere is approx- 
imately 0.1 km”’ (Reference 17). 

The validity of the exponential model for purposes of radio refraction correction has been well 
established (References 7, 12, and 18). 


With regard to Equation 1 5, the optical surface refractivity can be calculated from 




11.6P 

T 


(16) 


where 


P = total atmospheric pressure (millibars) 

and 

T= air temperature (degrees Kelvin). 

The next step is to determine what the Minitrack direction cosine measurements should be, based 
upon the line-of-sight angle to the aircraft (E + e of Figure 8 ). During the calibration, the refracted 
starlight is made coincident with the flashing light on the aircraft. Let cos a and cos |3 be the line-of- 
sight direction cosines to the star as obtained from a star catalog, and let cos a' and cos |3' be the de- 
sired line-of-sight direction cosines to the aircraft. Then 


cos a' = cos a 



( 17 ) 
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and 


From previous definitions. 



and 


cos a = sin AZ cos E , 
cos |3 = cos^Z cos E , 

cos a' = sin AZ cos (F + e) = sin AZ (cos E - e sinE) , 


cos ]3' = cos AZ cos (E + e) = cos AZ (cos F - e sin F) . 

Thus the ratios cos a'/cos a and cos jS'/cos j3 (Equations 17 and 18, respectively) are given by 

cos a' cos|3' , ^ ^ 

= = 1 - e tan E = (I - K) . 

cos a cos p 

By combining Equations 15 and 21, the correction factor (1 - becomes 


( 18 ) 


(19) 

( 20 ) 


( 21 ) 


and 


(1-70 = 


A^„(10"^)[l 


( 22 ) 


cos a' = cos a{\- K) 

(23) 

cos |3' = cos /3(1 - , 

(24) 


cos a and cos (3 being line-of-sight direction cosines to a star coincident with the flashing light and 
cos a and cos jS' line-of-sight direction cosines to the aircraft. (Equation 22 is currently employed in 
the calibration calculation at the Physical Science Laboratory.) 


The next section will show that the expected radio phase meter readings Aip^ and A(p^ are given by 


and 


= [b cos a(l - K)n] 
Av?^ = [b cos j3(l - K)n] 



(25) 

(26) 


where 


b = baseline length, 

cos a and cos /3 = direction cosines corresponding to line of sight to star occulted by the 
flashing light. 
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( 1 - ^ = correction factor to obtain line of sight to aircraft. 


and 



n(Ji)dh = average radio index of refraction between tracking 
station and aircraft altitudes, 


= free-space wavelength of received frequency. 


5.2 Calibration and Radio Refraction 

During the aircraft calibration flyby at a given site, radio phase data are recorded and the differ- 
ence between the observed phase data and that “expected” is sent to Goddard Space Flight Center as 
a correction to be applied to subsequent raw Minitrack phase data. During this calibration the expected 
Minitrack radio phase measurement is based upon the precise line of sight to the aircraft adjusted for 
prevailing radio refraction conditions. The previous section discussed the conversion of star line of 
sight to aircraft line of sight. This section will derive the radio refraction correction factor. 

During calibration the ionosphere plays no part since the aircraft is necessarily well below the 
ionosphere’s 85-km lower altitude limit. A most important point is that the radio phase readings 
measured during calibration are compared with precise line-of-sight data altered slightly to take into 
account the tropospheric radio refraction. The calibration is not a means for later adjusting raw space- 
craft tracking data to either angle of arrival or line-of-sight measurements; the calibration accounts for 
on-site biases associated with the Minitrack antenna system and electronics and does not provide a 
refraction correction for subsequent raw data recordings. 

In order to calculate what the interferometer should read for a given line-of-sight geometry, it is 
necessary to calculate the electrical path difference corresponding to the geometric path difference 
- iSj (Figure 9). It will be shown that, for the aircraft geometry, the corresponding electrical path 
length is given by - S.^), where n is the average radio index of refraction between tracking site and 
aircraft altitudes. 

With reference to Figure 9, in the absence of an atmosphere the free-space interferometer phase 
measurement would be given by 




27t(5'2 


(27) 


where 

iSj = geometric distance from the phase center of the antenna at one end of the baseline to the 
aircraft, 

S., = geometric distance from the phase center of the antenna at other end of the basehne to 
the aircraft, 

and 

X^ = free space wavelength = c//^, c being the speed of light in vacuum and the frequency of 

the received radio energy. 
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H» = APPROXIMATE AIRCRAFT HEIGHT (NEGLECTING 
EARTH CURVATURE) 

Hs = STATION ALTITUDE 
Figure 9— Aircraft calibration geometry. 


The frequency will differ from the transmitted frequency by the one-way Doppler shift (r/c), 
where r is the radial velocity of the aircraft beacon relative to the station. At 136 MHz, this is en- 
tirely negligible relative to the expected uncertainty in the measured surface index of refraction of 
+\N unit; that is, the ratio of ISN (±1 ) to (~300) is always large compared to rjc. The physics link- 
ing surface refractivity to wet- and dry-bulb temperature measurements and statements of expected 
accuracy is covered in detail in Reference 7. 

Again with reference to Figure 9, the actual interferometer phase readings are given by 



n{s)dS2 


(28) 


As pointed out by Bean and Thayer (Reference 14) of the National Bureau of Standards (NBS), 
Equation 28 can be integrated assuming flat-earth conditions for elevation angles > 10°, slant range 
r<a^ sin {a^ = earth radius), and baselines of a few hundred meters or less. 

These conditions are applicable to the aircraft flyby and hence, again with reference to Figure 9, 
Equation 28 can be approximated by 
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(29) 


X„Vsin72 


dh. 


which neglects ray curvature. Therefore, by the law of sines. 


and 


sin 7j sm 7 


sin 


72 sin 7 


Equation 29 becomes 


A(p = 


271(5’, - 


,- 5 i)/ 1 X f^A 


277 ( 52 - 5 .) 
dh = ; n . 


(30) 


With reference to Figure 9, 52 - 5p at aircraft slant ranges, is given by 

52 - 5j = h cos i// . 

For one Minitrack baseline, 1// = a'; for the other orthogonal Minitrack baseline, = jS'. In relation to 
the star direction cosines, a and |3, Equation 30 can be written as 

277 

Ap^ = [b cos a(l - XOn] — 


and 


^[b cos 13(1 -K)h]^, 


(31) 


where the factor (1-/0 was defined by Equation 22. 


In Equation 3 1 , Ap^ and Ap^ represent the expected Minitrack phase readings during the aircraft 
flyby when the flashing light appears coincident with a star whose line-of-sight direction cosines are 
cos a and cos (3 respectively. 

'The average radio index h is calculated using the NBS exponential profile n(h) = 1 + 

This profile, when integrated and averaged, results in 


n = 1 


A^^(10“^) [1 - 


(32) 


The decay factor k is calculated with an empirical relation formulated by NBS based upon a fit to 888 
sets of data (Reference 12) 
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k = log^ 






where AN represents the change in refractivity in the vertical direction 1 km above the tracking station 
(Reference 12): 




(33) 


Finally, the surface refractivity is calculated from surface meteorological measurements with the 
relation 




77.6 


P + 


4810C/RH)' 


(34) 


where 

= 1 + Nj(10”^), the surface radio index of refraction, 

T = air temperature (°K), 

P = total air pressure (millibars), 

Cj = saturation vapor pressure (millibars), 
and 

RH = relative humidity (Reference 1 9). 


The optical surface refractivity as previously indicated by Equation 16, is Equation 34 with- 
out the water vapor term. 


Equations 3 1 through 34 are used by New Mexico State University in their reduction of Mini- 
track aircraft calibration data for Goddard Space Flight Center. 


5.3 Angle of Arrival Correction Equations 


As a matter of interest, the foregoing principles can now be applied to derive the flat-earth ap- 
proximation for the angular difference between elevation angle of arrival and line of sight for both 
optical and radio refraction. The only difference between the optical and radio cases is that n for the 
radio refraction includes the effect of water vapor within the earth’s atmosphere. 


For an interferometer such as Minitrack, which has intersecting orthogonal baselines running east- 
west and north-south, from previous considerations, the phase readings during a spacecraft track are 
given by 


^ a 


2irb cos E _ 

yi 


sin AZ 


and 




2nb cos E _ 


n cos AZ . 


(35) 
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Here, E represents the line-of-sight elevation angle to the source and AZ the azimuth angle. The 
phase measurement is also represented by the phase delay in terms of the number of wavelengths of 
delay within the media surrounding the antenna. The effective wavelength within the media in the 
vicinity of the antenna is In terms of the elevation angle of arrival E^, the following equations 

can thus be written: 


and 


2-nb cos E^ 

Ap = — ; — sin AZ 

“ KM. 

O' S 


2-nb cos E^ 

Ap„ = ; ] cos AZ . 

P 

o' S 


Equations 35 and 36 combined yield 


(36) 


n cos E = cos E^ . (37) 

This equation links elevation angle of arrival E^ to line-of-sight angle E ; but E^ = E + EE, where AE is 
a small perturbation, and Equation 37 can be expanded to obtain the approximation 

A£= ^1 -^^cot£:. (38) 

Here AE is the elevation angle-of-arrival error due to the atmosphere. This result was also derived by 
R. H. Paul (Equation 19 of Reference 14) using a somewhat different approach. Equation 38 can also 
be written in terms of refractivity N since 

n= 1 -l-iV(10~^) 

and 

1 +7V^(10'6). 

In all cases the deviation from unity is 10“^ or less; thus. Equation 38 can be written as 

AE = iN^~ N) cot £■( 10"®), (39) 

where, as before, represents the radio surface refractivity varying from 200 to 400, depending upon 
geographical location and local meteorological conditions, and N is the average refractivity ^etween 
the station altitude and the altitude of the signal source. If the source is a spacecraft, then N includes 
the effect of the ionosphere, where the absolute sign of the refractivity term is negative (Reference 6). 
Equation 39 thus also indicates how the presence of the ionosphere increases the total ray bending as 
suggested by Figure 5. If the altitude of the spacecraft increases much above 5000 km, or for all _ 
spacecraft operating at 2 GHz or higher such that the ionosphere is no longer a significant influence, N 
tends to zero and Equation 39 reduces to 

AE = N^il0~^) cotE , 
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and, to a first order, 


cot E - cot E„ , 

o ’ 

/SE ^N^{\QT^) cot E^ . 


or 


( 40 ) 


This result, presented earlier as Equation 2, is the usual tropospheric-bending correction applied to 
elevation- angle spacecraft tracking data for elevation angles above 10°. The same result is used for 
optical correction to star-light angle of arrival, except that the optical surface refractivity, is used 
instead of N^, the radio-surface refractivity. 

Finally, Equations 38 and 40 can be used to derive Equation 14, presented as the angular differ- 
ence e between the line of sight to a star coincident with a flashing light on the Minitrack calibration 
aircraft and the line of sight to the aircraft. 

With reference to Figure 8, the angle e is the difference between the total star-light bending 
A^q( 10”^) cot E and the bending of the flashing-light ray path relative to the line of sight to the aircraft 
(1 - nln^) cot E: 


or 


(«„-!)- 



cot E , 



cot E . 


(41) 


Equation 41 was presented as Equation 14 and presents a good approximation of elevation angle 
distortion for elevation angle E^ greater than 1 0° and slant ranges to the spacecraft much less than 
sin E^, where is the earth radius. 


6.0 EFFECT OF IONOSPHERE 

ON MINITRACK MEASUREMENT 

The elevation angle bias present in Minitrack measurements below elevation angles of 45° can 
largely be attributed to residual error after nominal ionospheric corrections have been applied. Unlike 
the troposphere, the ionosphere is much less predictable, and, as mentioned previously, ionospheric 
refraction effects vary inversely as the square of the radio frequency. One of the major sources of 
ionospheric modeling error has to do with the fact that the radio signal often traverses portions of the 
ionosphere several thousand kilometers from the tracking site, a situation where the usual, spherically 
stratified Chapman model does not strictly apply. One means of minimizing the ionospheric refraction 
error is to use only data where the signal is traversing local night regions. Recent studies, however, in- 
dicate ionospheric propagation anomalies peaking at local midnight at the South American Minitrack 
sites (Santiago, Lima, and Quito) (Reference 20). These equatorial region anomalies appear to pre- 
dominate in the east-west baseline which, as stated in Reference 20, is consistent with existing theories 
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of ion plasma density alignment north- to-south along the earth’s magnetic field lines. Study of the phe- 
nomena and the effect on NASA 136-MHz telemetry, tracking, and command signals is continuing. *t 

6.1 The Interferometer Refraction Error Equation 

The inherent tropospheric correction to interferometer data was discussed in section 4.0. It was 
shown that the elevation angle calculation is given by 

£ -A(10"^)cot£ 
o ' o 

and (42) 

E' =E^-N^ilO~^) cot E , 

forE^ >10°, where, 

E^ = elevation angle of arrival, 

E' - interferometer calculated angle, 

and 

E = line-of-sight elevation angle. 

It was also shown in section 5.3 (Equation 38) that the line-of-sight elevation angle E is linked to 
angle of arrival elevation E^ by 


cot Ej . (43) 

Thus the total interferometer error is found by combining Equations 42 and 43 to obtain 

Interferometer Error ~ E' ~ E cot E , (44) 

which corresponds to Bean and Thayer’s Equation 26 in Reference 14. Equation 44 is a good approxi- 
mation for E> 10°, interferometer baselines of a few hundred meters or less (such as Minitrack), and 
slant range to spacecraft r sin E, where is the earth radius. Equation 44 is useful in qualita- 
tively describing the effect of the ionosphere on Minitrack data when the spacecraft is within the iono- 
sphere ( 1 00 km < altitude < 1 000 km). At 1 36 MHz the contribution to N due to the troposphere is 
negligible compared to the ionospheric contribution. This can be shown as follows. According to 
Equation 32, for altitudes significantly above the troposphere scale height (~7 km), the tropospheric 
average refractivity is given by 


E = E^ 


1 -- 


N. 


N 
^ 

kH' 


(45) 


*Golden, T. S., “A Proposal for Measuring the Correlation Distance of Kqiiatorial Ionospheric Disturbances,” NASA-GSI'C 
Document X-5 25-68-302, August 1968. 

tGolden, T. S., “Effects of Differential Signal Amplitudes in the Minitrack System,” NASA-GSl’C Document X-525-68-14, 
January 1968. 
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where 


H = spacecraft height above the tracking station 

and 

k = decay factor in exponential model, 1 /A: = 7 km. 

For/f = 200 km, = 12, which is small when compared to the ionospheric contribution. The 
ionospheric refractivity is given by (Reference 6) 

-AQ3N(h) 

Nih) = ;^(106), (46) 

/ 


where 


N - electron density (electrons/meter^) 

and 

/ = frequency (Hz). 

For a daytime peak density at the height of maximum ionization, takes on values of typically 
0.5 X 10^^ electrons/m^, and at 136 MHz this corresponds to a maximum ionospheric refractivity of 
the order of -2000. This value will taper off by a factor of 10 at the lower edge of the ionosphere 
(H = 85 km). However, the average contribution will be significantly greater than the tropospheric 
contribution = 1 2. An interesting observation is that, again by Equation 44, and noting that the 
average refractivity due to the ionosphere is negative, the ray bending will be in the same direction as 
that due to the troposphere. With the interferometer, of course, the tropospheric bending is already 
inherently corrected for by the usual angle calculation discussed in section 4.0. 

The flat-earth model (Equation 44) is useful for qualitative ionospheric refraction discussion only 
when the range of validity includes the constraint r<a^ sin E. However, a more satisfactory equation, 
which includes earth curvature, has been presented by Bean and Thayer, where the interferometer 
error is given by 


Interferometer Error = 


f 


lO'^cotE'/ N(h)dh 
’0 


r sin E sin^ E^ 


for E> \Q° and r > b{b < 50 km), where 
E = line-of-sight elevation angle, 

N{h) - vertical refractivity profile including both troposphere and ionosphere, 
r = slant range station to spacecraft, 

H = altitude of spacecraft above station, 


(47) 
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= earth radius, 

and 

b = interferometer baseline (6 = 120 meters for Minitrack). 

For the flat-earth approximation (j = H and sin E), Equation 47 reduces to ~N{ 10“^) 

cot E, presented previously as Equation 44. 

Equation 47 can be derived from Snell’s law for a spherically stratified media. The result can also 
be obtained by ray tracing through the atmosphere (again using Snell’s law) to determine the differ- 
ence between angle of arrival and line of sight. This can then be linked to the interferometer error by 
subtracting out the first order tropospheric correction of Af^(10~^) cot E. 

6.2 Minitrack Ionospheric Refraction 

Tliis section presents the Minitrack bias error that can be expected if ionospheric effects are ne- 
glected. A nominal correction based upon NBS monthly worldwide peak electron density predictions 
and some mathematical distribution model (such as a Chapman profile) often reduce the error by only 
50%. The greatest uncertainty is related to the F^-layer vertical-sounding (i.e., zero) maximum usable 
frequency (MUF). The peak electron density, and hence refractivity, is proportional to the square of 
the F 2 (zero) MUF. 

As indicated by Figure 5 and Equation 47, the Minitrack refraction error is separable into tropo- 
spheric and ionospheric components. At frequencies such as 136 MHz, where ionospheric bending is 
appreciable, the tropospheric contribution in Equation 47 is generally negligible if elevation angles 
above 10° and spacecraft altitudes above 100 km are considered. During aircraft calibration the re- 
verse is true and only the troposphere is a consideration. By ray tracing (Reference 21 ) through the 
ionosphere and troposphere separately and then comparing with a trace through ionosphere plus tro- 
posphere, one may verify the validity of superposition of bending. The difference between total bend- 
ing due to the troposphere alone and Aj(10“®) cot E^ comprises the tropospheric contribution to the 
error in Equation 47. The ionospheric-caused elevation error in Equation 47 can be obtained directly 
from a ray trace through the ionosphere alone. 

Table 1 is an example of the Minitrack refraction-error magnitude associated with a track of the 
2000-km-altitude geodetic satellite GEOS 1 which was launched on 6 November 1965. It should be 
emphasized that exact refraction values will vary from site to site and depend upon the time of day as 
well as time of year. For the table a peak ionospheric density of 0.8 X 10^^ electrons/m^ was assumed 
for daytime and 10^^ electrons/m^ for nighttime. 

In Table 1, the term bias refers to that value which must be subtracted from the measured eleva- 
tion to obtain true elevation. The minus sign on the troposphere bias simply indicates the extent to 
which the inherent A^( 10“^) cot E correction of the interferometer overcompensates for the tropo- 
spheric bending. Also, the ionospheric refraction is a maximum for satellites at an altitude of approxi- 
mately 500 km. For an altitude of 500 km, the ionospheric error values of Table 1 will be increased 
by approximately 25% at 10°, 50% at 20°, and 70% at 30°. The ionospheric bias tabulated in Table 1 
is also plotted in Figure 10. The parameters used to obtain the results of Table 1 and Figure 10 via a 
ray trace are as follows: 
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ELEVATION BIAS (MILLIRADIANS) 


Table 1— Minitrack refraction bias (no corrections applied; satellite altitude = 2000 km). 


Elevation angle 

Tropospheric bias 

Ionospheric bias 
(mrad) 

(deg) 

(mrad) 

day 

night 

10 

-0.06 

2.25 

0.30 

15 

-0.02 

1.65 

0.20 

20 

-0.01 

1.25 

0.15 

30 

-0.005 

0.80 

0.10 

40 

-0.002 

0.50 

0.05 

60 

-0.001 

0.25 

0.02 

80 

~ -0.0005 

0.10 

0.01 


MIIMITRACK IONOSPHERIC PARAMETERS 


(JANUARY 1966 PREDICTION) 



REFERENCE 22 



DAY 

NIGHT 

MAXIMUM ELECTRON 
DENSITY (ELECTRONS/M 3) 

0.8(I0'2) 

0.1 (I0'2) 

HEIGHT OF MAXIMUM (KM) 

300 

250 

SCALE HEIGHT (KM) 

83 

66 

SATELLITE HEIGHT 

-2000 KM 



ELEVATION ANGLE (DEGREES) 


Figure 10— Minitrack elevation angle of arrival bias due to ionosphere if no correction applied. 


TROPOSPHERE: 


NBS exponential model 

7Vj = 350 k = 0.16km"i 
n{h)^ 1 +M/t)(10"^) 

0</z<30km (48) 
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IONOSPHERE: 


Chapman Model 

-A03N(h) 

N{h) = ^ 10^ 

/= 136 MHz 

n(,h)= \+N{,h)\QT^ /z > 85 km 

" 7 V^el/ 2 (l-z-e-^) (Reference 23) (49) 

= maximum electron density = 0.8(10^^) electron/m^ representing nominal daytime value. 


H ’ 

= 300 km representing a typical daytime height of maximum ionization, 

= scale height == 83 km, 

= 2000 km (GEOS 1 altitude). 

For a more precise determination of refraction effects, all of the foregoing parameters must be 
estimated in view of particular station geographical locations, time of year, time of day, phase of sun 
spot cycle, prevailing meteorological conditions, and so on. It should be pointed out that at present, 
except for special tracking tests, only Minitrack data above elevation angles of 80° are used for orbit 
computation. 

Figures 1 1 and 12 are plots of elevation angle residuals (observed elevation minus orbit calculated 
elevation) when no refraction corrections were applied to a 7-day stretch of GEOS 1 Minitrack data 
(31 December 1965 to 6 January 1966). It is felt that the trended portion of these curves can be 
attributed primarily to refraction within the ionosphere. 

Finally, a most important point is that the tropospheric error indicated in Table 1 would be the 
only significant atmospheric refraction bias for an interferometer operating at frequencies above 2 
GHz. For such a situation, and for satellites where the slant range approaches the magnitude of the 
earth radius multiplied by sin E, the interferometer refraction bias is given by the term associated with 
1/a^ sin^ E of Equation 47. This is the expected interferometer refraction bias presented by the 
National Academy of Sciences (Reference 13) for interferometers with a baseline less than 3 km and 
operating at high enough . frequencies such that the ionosphere is of no concern. The results in Refer- 
ence 13 check those for the tropospheric refraction presented in Table 1, the latter having been 
obtained from ray trace considerations. 

One might wonder why 136 MHz is used instead of perhaps 2 GHz for the Minitrack frequency 
to eliminate ionospheric refraction. The primary reason is historical because in 1956, when Minitrack 


N.. 


and 


Z 

H 

h 
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Figure 1 1 —Minitrack elevation angle orbital residuals 
GEOS 1 (no ionospheric correction). 


ELEVATION RESIDUAL 
(MILLIRADIANS) 


1.5 


FORT MYERS (FLORIDA) 
l5H<t<l9H 
LOCAL TIME 




COLLEGE (ALASKA) 
I7H< t<20H 
LOCAL TIME 


Figure 1 2— Minitrack elevation angle orbital residuals 
GEOS 1 (no ionospheric correction). 
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was first implemented, efficient S-Band spacecraft-qualified hardware was not available. However, the 
use of 136 MHz stUl has certain advantages over S-Band, such as ease of signal acquisition and ambigu- 
ity resolution. That is, the transfer of radio energy from a spacecraft using a low gain antenna to a 
given ground-located effective antenna aperture is independent of operating frequency (Reference 24). 
However, the respective beam and lobe widths vary inversely as the operating frequency. The antenna 
patterns corresponding to the Minitrack fine resolution array* are given in Figure 13. Note that there 
are two such arrays for each of two orthogonal baselines. The complete Minitrack antenna layout as 
given in Reference 1 5 is shown in Figure 14. 

6.3 Ionosphere Model Parameters 

The Chapman model presented in section 6.2 as Equation 49 requires three input parameters: 
the maximum electron density N^, the height of the maximum h^, and the scale height H. Of these 
three, the most critical parameter is the maximum density. The scale height has been formulated in 
terms of (Reference 25); 


//= 1 . 66[30 + 0 . 2 (/ 2 ^- 200 )] , ( 50 ) 

where H and are in kilometers. 

The height of the maximum, if no other information is available, can be taken as h^ = 300 km 
for daytime and = 250 km for nighttime (Reference 26). The overall refraction is not a critical 
function of and H. 

The following indicates a means for estimating the peak electron density of the ionosphere, uti- 
lizing the Institute for Telecommunication Sciences-Environmental Science Services Administration 
(formerly Central Radio Propagation Laboratory-NBS) monthly publication. Ionospheric Predictions, 
in conjunction with a worldwide, gyro frequency chart. 

Ionospheric Predictions presents Fj (zero) MUF worldwide charts (Reference 22) for 3 months 
in advance. The charts present monthly mean values in 2-hr increments, since the ionization variations 
are primarily diurnal during any given month. The “zero” in F 2 (zero) refers to zero ground “skip 
distance,” i.e., vertical incidence ionosphere sounding. It can be shown (see for example. Reference 
27) that for vertical incidence the critical frequency (highest frequency signal that will be returned to 
earth) is a function only of peak electron density and magnitude of the earth’s magnetic field in the 
region of reflection. The direction of the earth’s magnetic field does not influence the vertical critical 
frequency. Because of the presence of the earth’s magnetic field, three discrete critical frequencies 
exist: the ordinary critical frequency (this frequency would be the only critical frequency in the ab- 
sence of the earth’s magnetic field) and two so-called extraordinary critical frequencies. The F 2 (zero) 
MUF charts are for the highest critical frequency, or the upper extraordinary frequency. This upper 
frequency results in the more predominant reflection since signals at the lower frequencies encounter 
appreciable attenuation. 


'“'Lantz, Paul A., and G. R. Tliibodcau, "NASA Spacc-Directcd Antennas," NASA-GSi'C Document X-525-67-430, September 

1967. 
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Figure 14— Relative location and function of antennas 
in a Minitrack antenna field. 


The peak electron density is obtained from the ordinary critical frequency as 


80.6 


( 51 ) 


where 

/ = ordinary critical frequency (Hz) 

and 

- peak electron density (electron/m^). 

The p 2 (zero) MUE chart is for the upper extraordinary frequency (e.g., Figure 15) and is related 
to the ordinary critical frequency as follows: 







( 52 ) 
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60° 90° 120° 150° 180° BO® 120° 90° 60° 30° 0° 30° 60° 

EAST WEST WEST EAST 



EAST WEST WEST EAST 

PREDICTED MEDIAN MUF (ZER0)F2 (MHz) 


(PLOT OF fy ) 

Figure 1 5— Example of F 2 (zero) MUF prediction (Reference 22). 


where 

/ = upper extraordinary critical frequency as obtained from F, (zero) MUF chart, 
f - ordinary critical frequency required to calculate maximum N (Equation 51), 

and 

= gyro frequency, which couples the effect of the earth’s magnetic field into the iono- 
sphere index of refraction (Figure 15). 


The last quantity, /^, is defined (Reference 27) as 


where B is the earth’s magnetic field (Gauss). Since /^/4 < in most cases, 
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SOUTH NORTH 




or, the ordinary critical frequency is simply 


fh 


■4 4ci 2 ■ 


(53) 


The predicted and observed Fj (zero) MUF values seldom differ by more than 1 or 2 MHz, with 
the larger differences occurring at the beginning or end of a given month. Unfortunately, the maxi- 
mum electron density varies as the square of the critical frequency (Equation 51); so for a value of 
/^ = 5 MHz, under worst case conditions (error of 2 MHz) the electron density could be off by a factor 
of 2. If the predicted critical frequency is off by only 1 MHz, then the maximum electron density in 
this case would be correct to within 30%. During mid-month, the F 2 (zero) MUF predicted values are 
accurate to within 80% to 90%. 


The maximum electron density is given by 




80.6 


(54) 


where 

f = frequency for appropriate geographical location, month, and time of day, as obtained 
■^1 

from an F 2 (zero) MUF chart such as Figure 15, 

4 = gyro frequency appropriate to a particular tracking site (Figure 16) (note that the gyro 
frequency is not a function of time), 
and 

Nj^ = estimated maximum electron density in the region above a particular tracking station 
(electrons/m ^). 

The maximum electron density calculated by Equation 53 is used in the Chapman model 
(Equation 49) as a scale factor for the vertical electron density profile. To obtain best results, the pro- 
file should be calculated for the region where the radio energy intercepts the region of maximum den- 
sity. Generally, however, only the region above the tracking site is considered and spherical symmetry 
is assumed. The error thus introduced depends on satellite altitude and the relative steepness of con- 
tours such as given by Figure 15. 

7.0 CONCLUSIONS 

The atmosphere, for purposes of radiowave propagation, can be separated into two parts: the 
troposphere (sea level to 30 km) and the ionosphere (85 km to 1000 km). Troposphere refraction 
effects are frequency independent in the range 100 MHz to 20,000 MHz. Ionospheric refraction is 
frequency dependent varying inversely as the square of the operating frequency. At 136 MHz iono- 
spheric refraction is significant. The atmosphere influences Minitrack measurements during aircraft 
calibration and during spacecraft tracking. 
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Figure 16-Gyro-frequency map of the world (Reference 28). 


AIRCRAFT CALIBRATION 

1. Only the troposphere is a consideration for both radio and optical (star background and flash- 
ing light) refraction. 

2. The NBS exponential atmosphere coupled with flat-earth geometry is satisfactory for model- 
ing aircraft calibration refraction. 

3. For this model only the scalar quantity, surface refractivity (optical and radio), is required. 

SPACECRAFT TRACKING 

1 . Only the ionosphere is a consideration since the troposphere is inherently corrected for by the 
Minitrack calculation. 

2. The ionosphere can be modeled using a Chapman-type profile to obtain estimates of iono- 
spheric induced bias. 
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3. Recent top-side satellite soundings suggest that electron density decays with altitude some- 
what less rapidly than predicted by the Chapman model. An exponential decay in electron 
density above the height of maximum density appears to approximate more closely the true 
profile. 

4. The parameters maximum electron density, height of maximum density, and scale height are 
required. 

5. The extent to which corrections are successful is highly dependent on validity of estimates for 
maximum electron density. 

6. Published prediction for maximum density used at the beginning and end of a given month 
can be such that as much as 50% of ionosphere residual bias remains after corrections are 
applied. 
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